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lO [ Abstract 

We study the gravitational, electromagnetic and scalar field perturbations on the near-horizon 
geometries of the even-dimensional extremal Myers-Perry black holes. By dimensional reduction, 
the perturbation equations are reduced to effective equations of motion in AdS2- We find that some 
modes in the gravitational perturbations violate the Breitenlohner-Freedman bound in AdS2- This 
result suggests that the even-dimensional (near-)extremal Myers-Perry black holes are unstable 
against gravitational perturbations. We also discuss implications of our results to the Kerr-CFT 
correspondence. 
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00 

§■ 1 Introduction 

00 

q " Black hole solutions have been fundamental objects in the study of general relativity theory. Higher- 
■ dimensional spacetimes introduced in some particle physics models motivated to consider their gen- 
eralizations to higher dimensions, and as a result diverse solutions of black objects were found [TJ [2]. 
A feature specific to those higher-dimensional black hole solutions is that they may be unstable, even 
^ I though the Schwarzschild and Kerr black holes in four dimensions are stable. Those instabilities signal 
c$ ' appearance of new blanches of black hole solutions, and also tell us that some other solutions will be 
realized as a result of time evolution. As such, the stability analysis of the higher- dimensional black 
objects is one of the most important steps to understand the phase structure of the solutions in the 
higher dimensions and also the final states of the higher-dimensional spacetimes. 

The generalizations of the (four-dimensional) Kerr black holes in d- dimensional spacetime are the 
Myers-Perry (MP) black holes [3J, which may have \_(D — 1)/2J different angular momenta. By an 
intuitive argument, it was predicted that "ultraspinning" MP black holes should suffer from an insta- 
bility [3]. Here, by ultraspinning, we mean that Hessian matrix, H a p = —d 2 S/dx a dx b (x a = M, Jj), 
for a black hole with mass M, angular momenta J, and horizon entropy S has at least two negative 
eigenvalues [5] Q Refs. O d El E] have clarified the conditions for the ultraspinning instability of the 
singly-spinning MP black holes, and also Refs. [3 EDI Ell E21 E3] studied the odd-dimensional MP 



1 In Ref. [5], it was shown that any vacuum stationary black holes have at least one negative mode in the Hesse 
matrix. It is a manifestation of the property that black holes have negative specific heat. Ultraspinning MP black holes 
have at least one more negative mode in addition to that. 
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black holes with equal angular momenta or all but one angular momenta are equal. Since also the 
even- dimensional MP black holes with equal angular momenta have ultra-spinning regime [5], we can 
expect the instability to occur even for these black holes. The instability of these even- dimensional MP 
black holes is, however, yet to be shown so far. A part of the reason is that the geometries of the even- 
dimensional ones are cohomogeneity-2 and less symmetric compared to those of the odd-dimensional 
ones, which are cohomogeneity-1, and then the stability analysis becomes more involved. 

One of our purposes is to fill this gap in our knowledge about the MP black holes by assessing 
the stability of the even- dimensional MP black holes with equal angular momenta. Since the stability 
analysis of the full spacetime requires to solve complicated partial differential equations, it is desirable 
that we have more simpler ways to obtain implications of the instability. Ref. [14] initiated a study 
on a new analysis method for such a demand. Their method is based on a generalization of the Geroch- 
Held- Penrose formalism to higher dimensions pj)]. They chose some combinations of the Weyl tensor 
components as perturbation variables, and found that the variables satisfy decoupled equations of mo- 
tion when the near-horizon geometries of extreme black holes are taken as the background spacetimes. 
Their method can be regarded as higher- dimensional generalization of the Teukolsky formalism [16] , 
although it is applicable only to the Kundt spacetimes|§ In Ref. [T2j, it was conjectured that the in- 
stability of the near-horizon geometries imply the instability of the original full spacetime. A proof for 
this statement is provided for the scalar field perturbations, and some pieces of evidence were found for 
the gravitational perturbations, especially on the background of the odd- dimensional MP black holes. 
In this paper, we assume this conjecture to be correct and study the stability of the even- dimensional 
MP black holes with equal angular momenta by examining the perturbations on their near-horizon 
geometries. 

We also have a motivation from the Kerr-CFT correspondence [17], which is a conjecture that 
near-extreme black holes are described by 2d CFTs. From the fall-off behavior of the perturbations on 
near-horizon geometries, we can read off the conformal weights of right sectors in the dual CFTs [121 
[T8l [T9| [20] . In four and five dimensions, some interesting properties of the conformal weights have 
been found for the operators corresponding to the axisymmetric perturbations. One of them is the 
integerness of the conformal weights: all the conformal weights for vacuum near-horizon geometries 
take integral values. Another one is the universality of the conformal weights: all the U(l) x U(l) 
symmetric five- dimensional vacuum near-horizon geometries with vanishing cosmological constant share 
the same sequences of conformal weig htsl This result is surprising because these five- dimensional near- 
horizon geometries contain dimensionless parameters on which the conformal weights may depend in 
principle. In addition to that, five- dimensional extreme black holes can have S 3 or S 2 x S 1 horizon 
topology [21] [23]. Nevertheless, those various near- horizon geometries share the common sequences of 
the conformal weights. This result suggests that there is a "universal sector" present in all the CFTs 
dual to extreme rotating black holes in five dimensions. These features have been observed only in four 
and five dimensions. It was also observed that the effective mass of the perturbative fields becomes 
non-integer but rational numbers for any parameters when the background near-horizon geometry is 
unstable [12] . We will examine if these properties persist even for the even- dimensional MP black holes 
we are focusing on. 

In our analysis method, which was used also in Ref. [12] , the perturbation equations of the near- 
horizon geometries are Kaluza-Klein reduced to an equation of motion of a massive charged scalar 

2 The Kundt spacetime is the spacetime admitting a null geodesic congruence with vanishing expansion, rotation and 
shear. All known near-horizon geometries are the Kundt spacetimes. 

In four dimensions, the near-horizon extreme Kerr geometry is the unique U(l) symmetric vacuum near-horizon 
geometry |211 122] and does not have any dimensionless parameters. Thus, the universality is trivial in four dimensions. 
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field on AdS2 spacetime whose effective mass is determined by the eigenvalues of the angular part 
of the perturbation equations. For the background of the even- dimensional MP black holes, those 
angular part equations are reduced to one-dimensional ordinary differential equations (ODEs) which 
are coupled to each other if we classify the perturbations into scalar /vector /tensor modes according to 
their transformation properties on CP n_1 contained in the induced geometry on the horizon. We solve 
those equations using the numerical technique employed in Ref. [18], which analyzed the conformal 
weights for the general near-horizon geometries in five dimensions. 

This paper is organized as follows. After introducing the near-horizon geometries and the perturba- 
tion equations in Sec. [21 we show the results of the stability analysis for the gravitational perturbations 
in Sec. [3j As a result, we find certain modes in these perturbations become unstable. In Sec. HJ we 
show the results for the scalar and electromagnetic field perturbations. These fields are stable in the 
full geometries, and we observe that the analysis of the near-horizon geometries gives results consistent 
with this fact. As a byproduct of our analysis, we find that the spectrum of the mass of the effective 
AdS2 field, which is given by the eigenvalues of the angular part of the perturbations, show a well- 
organized structure. We give comments on this property in Sec. [5j We also discuss the implications of 
our results to the Kerr-CFT correspondence in Sec. El especially about the integerness and universality 
of the conformal weights mentioned above. Sec. [7| is devoted to the summary and discussions, and 
appendices show the technical details of our analysis. 



2 Perturbations of near-horizon geometries 

We summarize the near-horizon geometries of the even-dimensional MP black holes with equal 
angular momenta, and also the perturbation equations on this background in this section. We also 
introduce the instability condition for the near- horizon geometries which is argued in Ref. [12] . 



2.1 Near-horizon geometry 

In our study, we focus on the extremal MP black holes in d = 2n + 2 (n > 2) dimensions with all 
the angular momenta set equal. Their near-horizon geometries can be written as a fibration over AdS2 
given by 



ds 2 = L{6) 2 (^-r 2 dt 2 + ^Pj + ® 2 {6)d6 2 + ^ 2 (6)g a(S dx a dx^ + $ 2 (6) (dfi + A- ^rdt) 2 , (1) 
where < 6 < ir, 

L 2 (9) = * 2 fn(0) e2(g) = a 2 /nW ^ ) = ^ & {g) _ 4nVsin 2 £ 



(2n-l) 2 ' w 2n-l' w 2n - 1 ' v; (2n-l)/ B (0)' 

1 ^ ' 

f n (9) = 1 + (2n - 1) cos 2 9, h* = == = n, 

ny In — 1 

<p is 27r-periodic and g a ^dx a dx^ (a, j5 = 1, . . . , 2n — 2) is the Fubini-Study metric on CP n_1 with Kahler 
form J = ^dA normalized to have R a p = 2ng a p. From the near-horizon metric, we can read off the 
induced metric on the horizon (the surface of constant t and r) as 

ds 2 = g^dx"dx v = Q 2 {6)d6 2 + $> 2 (6)g aP dx a dx p + $ 2 (#) (# + A) 2 . (3) 
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2.2 Perturbation equations 



Ref. [12] gave a prescription to obtain separated perturbation equations for general near-horizon 
geometries. We explain it briefly below. 

To describe the perturbations of the near-horizon geometry, we introduce a null basis {£, n, rrii}, 
where i = L(6){-rdt + dr/r)/y/2, n = L(6){rdt + dr/r)/y/2 and rrii are orthonormal spacelike vectors 
orthogonal to both I and n. We define Qij = C, 



' abed 



7a m b i i c mj and (fi 



F ah t a m\, where C abcd is the Weyl 
tensor and F a b is the electromagnetic field strength, and use them as the perturbation variables. As 
for the scalar field satisfying (□ — M 2 ~)ijj = 0, we use the scalar field itself as the perturbation variable. 
In this paper, we focus only on axisymmetric modes along direction: d^Qij = d^ipi = d^ip = 0. For 
these perturbations, we can take the separation ansatz as 



^ = X o(t, r)Y{y), ^ = Re [ X i{t, r)Y 4 (y)] , % = Re [x 2 (i, r)Y k 



(4) 



Then, the perturbation equation is separated as follows. The radial equations are given by massive 
charged Klein-Gordon equations in AdS2 with homogeneous electric fields 



[(V 2 - iq.A 2 



A s ]xs = 0, (5 = 0,1,2) 



(5) 



where V 2 is the covariant derivative on AdS2 {ds 2 = —r 2 dt 2 + ^) , A 2 = —rdt is the effective gauge 
field which appeared as a result of the dimensional reduction, and q s = is are the effective U(l) charges. 
The angular equations are given by 



O^Y = \ Y, (O^Y) t 

where the operators are defined as 
-V„ (L 2 V M F) + L 2 M 2 Y, 



(0™Y) 



A 2 K 



2 1 pvi 



O (0)y 

(OCDy) 



(6) 



(7) 



~V (L A V p Y fl ) + (2- -^KkA Y p + L 2 (r^ + l -Rg^ 



Y v 



+ (4 + 2 ( k - 4^))[mV,] - ^ (dL%k u ^ Y u , 



(° i2)Y ),u = -jf?" (^VpV) + (e- ^ 2 k P k P ) V + 2£ 2 {R MP + Rg Wp ) Y' W) - 2L 2 R^Y p(T 



+ 



(dk) Mp - T2 {d(L 2 ) A k) Mp + 2(k- d(L 2 )) {pl V p - 2 (k - d(L 2 )) p V 



(/-I 



Y p \ 



(9) 



where the covariant derivative V M and the curvature tensors such as R pvpa are defined with respect to 



the induced metric on the horizon (j3]). 

The solution of Eq. (jSJ) at large r behaves as Xs 



~ A± , where 



(10) 



Following the arguments of Ref. p2] , we shall call the near-horizon geometry to be unstable when A s 
violates the effective Breitenlohner-Freedman (BF) bound, i.e., \ s < —1/4, and A± becomes complex. 
Ref. [12] conjectured that the sufficient condition for the full black hole geometry to be unstable is 
that the near-horizon geometry is unstable against axisymmetric perturbations. This conjecture is 
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proved in the case of the scalar field, and it was suggested that a similar argument applies also to the 
gravitational perturbations. 

In this paper, we argue the stability of the even-dimensional MP black holes with equal angular 
momenta based on this conjecture. We solve Eqs. fl6]) analytically if possible and numerically otherwise 
to obtain A s . Applying the conjecture to the resultant A s , we will argue the stability of the full black 
hole geometries. 



3 Instability against gravitational perturbations 

We start our discussions from the gravitational perturbations on the near-horizon geometries. The 
gravitational perturbations can be decomposed by tensor, vector and scalar harmonics on the base space 
CP™" 1 , which are labeled by a principal quantum number k = 0,1,2, ... . By the mode decompositions, 
the eigenvalue equations (Q reduce to coupled ODEs which depend only on the 9 coordinate. Solving 
the ODEs, we obtain the eigenvalues for the operator 0^ 2 \ The computation is straightforward but 
too technical to show the details here. Hence, we describe only some important results in this section 
and we defer the details of the calculations to the appendices. 



3.1 Gravitational tensor modes 

We start from the simplest components, which are the gravitational tensor modes. The tensor 
modes are decomposed by tensor harmonics in CP n ~ , and as a result the eigenvalue equations ([9]) 
reduce to an ODE for a single coordinate 9 given by Eq. (|64p . For n = 2, the tensor harmonics do not 
exist on CP 1 = S 2 . For n > 3, we obtain an analytical expression for the eigenvalues of as 

x ^ 2 {K + „(, + !) 1} { + „(, + 1)} _ 2(,,' - 3, + 1) _ 2, (< = Mi2j ... ) (n) 

n(2n — 1) n(2n — 1) n 

where k = 0, 1, 2, ... is the principal quantum number which labels the tensor harmonics. The parame- 
ter a = =Fl separates the tensor harmonics into Hermitian and anti-Hermitian ones. See App. IB.3~Tl for 
the further details of the tensor harmonics. Another integer £ is the quantum number along 9 direction, 
which parametrizes the 9 dependence of Y^ u . This A2 is always non- negative, and does not violate the 
BF bound. Thus, there is no indication of instability in the tensor modes. 



3.2 Gravitational vector modes 

Next, we study the gravitational vector modes. The vector modes are decomposed by vector har- 
monics in CP n_ and their derivatives. For n = 2, the vector harmonics do not exist on CP = S 2 . For 
n > 3, we find the eigenvalues are given by 

2 {* + ,,(<+ !)}{« + n(<+ l) + 2} _ C 

n [2n — 1) 2n — 1 

where k = 0,1,2, ...,£ = £ Q , £ Q + 1, .. . and (C, £0) are the set of integers given by 

(C,4) = (4-2n,0), (4,-1), (4,0), (4, +1) . (13) 

The gravitational vector modes are described by four free variables, as we see in App. IB. 3. 21 This is 
why we have four kinds of the eigenvalues described by different integer sets (C, £$) as listed in Eq. ( fl~3|) . 
The eigenvalues are all positive, and hence no instability is implied. 
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3.3 Gravitational scalar modes 



Finally, we consider the gravitational scalar modes. We can expand the scalar modes by scalar 
harmonics in CP n_1 and their derivatives. As a result of mode decomposition, we obtain the coupled 
ODEs of ten variables in general cases (k > 1 and n > 2) as the eigenvalue equations to solve. In the 
cases of k = 0, 1 and/or n = 2, some of the unknown variables drop out and we have to treat those 
cases separately. In any case, we find that the eigenvalues are written in an unified expression as 

2 {K + „« + l)-l }{K + „« + l) } _C_ 

n(2n-l) V ; 2n-V K 1 

where k = 0, 1, 2, . . . and £ = £q, £o + 1, • • • • For a general case with k > 1 and n > 2, we have ten 
eigenvalues described by the integer sets (C, £q) given by 

(C,4) = (-2(n-l),0), (0,-1), (0,0), (0,+l), 

(2n,-2), (2n,-l), (2n,0), (2n,0), (2n,+l), (2n,+2). (15) 

Note that the constant C takes only three different values, that is, C = —2(n — 1), and 2n, and each 
of them is associated with one, three and six eigenvalues labeled by different Z , respectively. For the 
special values of k and n mentioned above, the eigenvalues are described by 

k = 0: (C,£ ) = (-2(n-l), +2), (0,+2), (2n, +2) (16) 

K = l )n = 2: (C7,4) = (-2(n-l),0), (0,0), (0,0), (0,+l), (2n,0), (2n,+l), (2n,+2) (17) 
« = l,n>2: (C7,4) = (-2(n-l),0), (0,0), (0,0), (0,+l), 

(2n,0), (2n,0), (2n,+l), (2n,+2) (18) 
«>l,n = 2: (C7,4) = (-2(n-l),0), (0,-1), (0,0), (0,+l), 

(2n,-2), (2n,-l), (2n,0), (2n,+l), (2n,+2). (19) 

Comparing Eqs. (TT6|) — (TT9|) to Eq. ( TT5|) . we notice that some of the eigenvalue sequences described by 
the integers listed in Eq. ( fl5l) drop out in the special cases with low k and n. We notice also that Iq 
increases in some of the remaining sequences, which implies that low £ modes existing in general cases 
are truncated in the special cases. 

From the expression of the eigenvalues, we can find the parameter region for the instability, that 
is, A2 < —1/4. The eigenvalues other than (C,£q) = (2n, — 2) in k > 1 and n > 2 are shown to 
be non-negative, and thus they do not imply instability. Now, let us examine the eigenvalues for 
(C,£q) = (2n, —2). For £ > —1, we can show that the eigenvalues are non-negative again. However, for 
the lowest mode £ = —2, we have 

_ 2(K-l)(«-2n) 
Aa " „(2n-l) • (20) 

This eigenvalue can be smaller than —1/4 depending on k and n. The instability condition can be 
written as 

K _< K < K+ , K± ^ n + 2 ± V2(2n-l)(3n-2). (21) 

We show the values of k± for some n > 2 (d = 2n + 2 > 6) in Table [TJ We find that a several numbers 
of k satisfy the instability condition for any n > 2. We find also that the parameter region for the 
instability widens as we increase the number of dimensions n. It would be interesting to search for 
instability in the full spacetime geometry corresponding to these unstable modes on the near-horizon 
geometries. 
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d 


6 


8 


10 


12 


14 


16 


18 


20 




1.28 


1.41 


1.54 


1.68 


1.81 


1.94 


2.08 


2.21 




3.72 


5.59 


7.46 


9.32 


11.19 


13.06 


14.92 


16.79 



Table 1: k± against the number of dimensions d. The near-horizon geometries are unstable against 
gravitational scalar perturbation satisfying k_ < k < k + . 



4 Scalar and electromagnetic field perturbations 

Having finished the study on the gravitational perturbations, we move on to the scalar and electro- 
magnetic field perturbations. These perturbations are guaranteed to be stable on the full geometry by 
the following argument. These fields satisfy dominant energy conditions and, in addition to that, there 
is a global timelike Killing vector in the full geometry. Then, we can construct an energy integral whose 
integrand is non-negative everywhere [251 126]. Hence, any instability cannot occur. Therefore, if we 
show the stability of the near-horizon geometry against these perturbations, we can give a non-trivial 
check of the conjecture in Ref. [12j . 

As for the scalar field perturbations, we may take the scalar field as the perturbation variable, and 
it is decomposed by the scalar harmonics on CP n_1 . We show the result for this case in Sec. 14.11 As for 
the electromagnetic field perturbations, the perturbation variables are decomposed into the vector and 
scalar modes on CW n ~ . We show the results for each mode in Sec. 14.21 and Sec. 14.31 respectively. We 
will see that any perturbations are stable, which is consistent with the conjecture mentioned above. 



4.1 Scalar field perturbations 

The scalar field can be decomposed by scalar harmonics in CP" - . As a result, we obtain a single 
ODE as the perturbation equation, which is given by Eq. ( 148|) . When the scalar field is massless 
(M = 0), we analytically find the eigenvalues to be given by 

x 2{« + n(l + l)-l}{« + n(l+l)} 2(w-l) 

A °- n{2n - 1) (22) 

where k — 0, 1, 2, . . . and i = 0, 1, 2 . . . . This Ao is non-negative for any k, i and n > 2, and thus the 
AdS2 BF bound (Ao > —1/4) is not violated. This is consistent with the stability of the scalar field 
perturbation on the full geometry. We found evidence that this property persists even for a massive 
scalar field (M > 0). See Sec. IB. II for details. 



4.2 Electromagnetic vector modes 

Electromagnetic fields can be decomposed by vector and scalar harmonics in CP n_1 . Here, we 
consider the vector modes. For n — 2, there is no vector harmonics on CP 1 . For n > 2, we analytically 
find the eigenvalues for the vector modes are given by 

2{ K + n « + l)}{ K + n « + l) + 2}_ 4-j, 
n(2n — 1) 2n — 1 

where £ = 0, 1, 2, . . . . These eigenvalues are always positive, and thus these modes are expected to be 
stable. 
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4.3 Electromagnetic scalar modes 



Finally, we consider the electromagnetic scalar modes, which can be decomposed by scalar harmonics 
on CP n_1 and their derivatives. As a result of the mode decomposition, we obtain coupled ODEs for 
four variables in a general case (k, > 1) as the eigenvalue equations to solve. When k = 0, the scalar 
harmonics becomes a constant on CP" -1 and the scalar-derived vectors vanish. Then, we have only 
two variables to solve. For both cases of k > 1 and k = 0, the eigenvalues are written by an unified 
expression as 

Xl ~ M2^T) ^ + 1) + 2^i' (24) 

where £ — £o, £q + 1, . . . . For k > 0, the integer sets (C, £$) given by 

(C,4) = (-2(n-l),0), (0,-1), (0,0), (0,+l). (25) 
For k > 0, they reduce to 

(C,£ ) = (-2(n-l),+l), (0,+l). (26) 

The eigenvalues for these parameter sets are non-negative, and instability is not implied. This result 
along with that in the previous section are consistent with the stability of the full geometry against 
the electromagnetic perturbations. 



5 Comments on the effective mass spectrum 

In this section, we give brief comments on the spectrum of the effective mass of the radial field on 
AdS2, namely, the relationships among the eigenvalues we have clarified in the previous sections. 

Comparing Eq. ( Tl~3|) to Eq. fl23|) . we notice that the eigenvalues of the gravitational vector modes 
corresponding to the first one in Eq. ( 1T3]) . (C, £ ) = (4 — 2n,0), coincide exactly with those of the 
electromagnetic vector modes given by Eq. fl23|) . Since the gravitational vector modes have also the 
eigenvalues with C = 4, the eigenvalue sets of the electromagnetic vector modes is completely covered 
by those for the gravitational vector modes. 

We also find similar relationships among the eigenvalues of the scalar perturbations of the various 
fields as follows. 

• Scalar field and electromagnetic scalar modes: 

Comparing Eq. (I24p with Eq. ( I22p . we find that the eigenvalues of the electromagnetic scalar 
modes coincide with the those of the scalar field perturbations if we set C = —2(n — 1). The 
first ones in Eqs. ( 125]) and fl26l) for the electromagnetic scalar modes share such a property, while 
the eigenvalue for (C,£) = (— 2(n — 1),0) do not appear in the electromagnetic scalar modes 
with k = 0. It means that the eigenvalues of the scalar field and electromagnetic scalar mode 
perturbations for C = —2{n — 1) are almost identical to each other for any k and n, and the only 
exception is the eigenvalues for £ = of the electromagnetic scalar modes with k = 0. 

• Electromagnetic and gravitational scalar modes: 

The eigenvalues corresponding to the first four integer sets in Eq. ( II 5p for the gravitational scalar 
modes, which have C = and C = —2(n — 1), coincide exactly with those of the electromagnetic 
scalar modes given by Eqs. (124"]) and (T25]) . When k — 0, 1 and/or n = 2, we find some of the 
eigenvalues of the electromagnetic scalar modes, especially those for low £ modes, are missing in 
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the eigenvalues with C = and C = — 2{n — 1) of the gravitational scalar modes, as we can see 
from Eqs. (jIBHP) - 

To summarize, we find the inclusion relations among the eigenvalues for the vector perturbations 
expressed as 

{Ai of EM vector with C = 4 - 2n} = {A 2 of grav. vector with C = 4 - 2n} , (27) 

{Ax of EM vector} C {A 2 of grav. vector} . (28) 

Similarly, if we do not care about a small numbers of the exceptions for small k and n mentioned above, 
we have the relations among the eigenvalues for the scalar perturbations given by 

{A with C = -2{n - 1)} = {A x of EM scalar with C = -2{n - 1)} 

= {A 2 of grav. scalar with C = —2{n — 1)} , 

{Ai of EM scalar with C = 0} = {A 2 of grav. scalar with C = 0} , (30) 

{Ao} C {Ai of EM scalar} C {A 2 of grav. scalar} . (31) 

If we care about the missing eigenvalues for k — 0, 1 and/or n = 2, the equality in Eqs. (129]) and (130]) 
should be changed into "d", and Eq. (13T]) needs to be modified accordingly. 

We also notice that the set of the eigenvalues for fixed C have simple organizations which are 
described by £q being integers around zero. This property is most noticeable in the eigenvalue list (fT5]) 
for the gravitational scalar modes. There is only one eigenvalue with C = — 2(n — 1), which is associated 
with £$ = 0. As for the eigenvalues which share C = 0, there are three of them, and they have £$ = — 1, 
and +1. Similarly, there are six eigenvalues which share C = 2n, and they have £ = —2, —1, 0, 
0, +1 and +2. The eigenvalue sets for the scalar field and electromagnetic scalar modes inherit this 
structure of the eigenvalues, and we can see a similar relationship to hold also between gravitational 
and electromagnetic vector modes. 

Similar relations for the eigenvalues can be seen also in the results of Ref. [12] [18]. The eigenvalues 
for various near-horizon geometries in five dimensions are studied in Ref. [18], and it was observed 
that the eigenvalues of the scalar field perturbations are included in those for the electromagnetic 
perturbations, and those for the gravitational perturbations cover both of them. Ref. [12] studied the 
eigenvalues for the odd-dimensional MP black holes with equal angular momenta. By rewriting their 
results in the form similar to Eq. ( Fl4]) . we find the the eigenvalues of the vector perturbations in that 
case is expressed by a unified formula 

x _ 2{ K + n(C' + i)}{K + n(C' + i)+2} | C 

n(n — 1) 2{n — 1) ' 



where n > 2 in this formula corresponds to iV + 1 in Ref. [12] . The eigenvalues of the electromagnetic 
vector modes are expressed by (C, C) = (— (n — 4), 0), while the eigenvalues for the gravitational vector 
modes are expressed by 

(C,C") = (-(n-4),0), ( 4 >-£)> ( 4 ' + £)- (33) 

The inclusion relation between those two modes and also the simple composition of the eigenvalues for 
C = 4 resemble to those in our case to some extent, that is, the eigenvalues for the electromagnetic 
vector modes appears as a part of the eigenvalues for the gravitational vector modes, and the eigenvalues 
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with C = 4 in the latter are described by some numbers C around zero which defer by one from each 
other. Similarly, the eigenvalues for the scalar field, electromagnetic and gravitational scalar modes are 
expressed by a single formula 



A 



2{K + n(C' + \) - l}{ K + n (C + \)} 
n(n — 1) 



+ 




C 



(34) 



where the constants (C, C) are given by 
Scalar field: (C, C) 



( 
( 



(n-2),0) 



(35) 



Grav. scalar 



EM scalar 



(C,C) 



( 
( 



(n-2),0), (o.-0' (°' + £) 
(n-2),0), (o,~), 



(36) 



(n-2),-l), (-(n-2),+l), (3n-2,0). 



(37) 



It is fair to say that the inclusion relations for the eigenvalues with C — — (n — 2) and 0, and also the 
composition of the eigenvalues with C = are similar to those in our cases. The only exceptions are 
the eigenvalues which appears only for the gravitational scalar modes, the last three sets in Eq. ( [37]) . 
They do not share the same constant C, while the counterparts in our cases (the last six eigenvalues 
in Eq. ([15]) with C = 2) did. Instead of that, the first, fourth and fifth eigenvalues in Eq. ( [37]) shares 
the same C = — (n — 2) in the case of the odd- dimensional MP black holes. 

The origin for the composition of the eigenvalues described above is unclear so far, though their 
simplicity tempts us to, naively thinking, suspect the existence of hidden background mechanism to 
generate it. It would be interesting to study the mathematical and physical origins of these properties 
studying, e.g., the perturbations of other near- horizon geometries. 

6 Conformal weight in the Kerr/CFT correspondence 

In any four- and five- dimensional vacuum near-horizon geometry, operators dual to the gravitational, 
electromagnetic and massless scalar field perturbations preserving rotational symmetry have integer 
conformal weights [12] [18] [19] [20]. Let us study if such a property persists in the case of the even- 
dimensional MP black holes we have studied. 

As we have seen in Sec. E] the near-horizon geometries have unstable modes in the gravitational 
scalar perturbations. Conformal weights become complex for the unstable modes, as we can see from 
Eq. ([TO]) . From the formulae for the eigenvalues we have shown in Sees. [3] and H] we find that the 
conformal weights for the most of stable modes take irrational numbers. For example, a stable grav- 
itational scalar modes for k = 0, n = 2 and (C,£q) = (2n, 2) have a sequence of conformal weights 
given by Hr — 1/2 = v / 201/6, \/33/2, \/417/6, .... This result suggests the following scenarios about 
the integerness of the conformal weights: the conformal weights become integers only in four and five 
dimensions, or they may take integer values even in higher dimensions while they do so only when the 
background black holes are stable. 

Below, we discuss another property of the conformal weights, the universality, mentioned in Sec. [U 
In any five-dimensional vacuum near-horizon geometry, operators dual to axisymmetric perturbations 
have universal conformal weights. That is, they do not depend on parameters nor horizon topologies 
of extreme black holes. We consider six dimensions for simplicity and check if the universality holds in 
this case. 
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To compare with the six-dimensional MP black holes, we consider the MP black string solution 
whose metric is given by ds 2 = ds 2 (MP^) + dz 2 where <is 2 (MP5) is the five- dimensional MP metric. 
By the dimensional reduction along ^-direction, the gravitational perturbation on the six-dimensional 
MP black string solution is decomposed into five- dimensional massless scalar, electromagnetic and 
gravitational perturbations on the background of ds 2 (MP 5 ). The conformal weights for those fields 
have been studied in Refs. Q2] , and it was found that all the eigenvalues for these perturbations 
are integers. On the other hand, the gravitational perturbations on the near-horizon geometries of 
six-dimensional MP black holes has rational but non-integer eigenvalues as we have seen in Sec. [31 This 
result suggests again that the universality is a property exists only in four and five dimensions or only 
for stable black holes. 

7 Summary and discussion 

We studied the gravitational, electromagnetic and scalar field perturbations on the near-horizon 
geometries of the even-dimensional extreme MP black holes with equal angular momenta. As a result 
of our study, we find that some modes in the gravitational scalar modes become unstable while all 
the others do not. This result implies that the instability for the corresponding modes in the full 
background geometries, assuming the claim of Ref. [12] to be correct. Ref. [13] raised a question if 
there is any (near-)extremal Myers-Perry black hole which is stable in d > 6, and our results suggest 
that there is no such in even dimensions if the angular momenta are set equal. If the black holes are 
unstable in the extreme limit, it is reasonable to consider that they are also unstable sufficiently near 
extremality. It is important to confirm this expectation analyzing the stability of the full geometry. 

Recently, a new method for demonstrating dynamical instability was established [271 [28] . They 
found an inequality whose violation implies black hole instability. The method based on this inequality 
makes the stability analysis dramatically easier, because the inequality can be evaluated only from 
initial data that describes a small perturbation of the black hole and we do not need to solve time 
evolution from that data. This method is a hopeful approach for demonstrating the instability of the 
full geometries of the even- dimensional MP black holes with equal angular momenta. 

For small angular momenta, the MP black holes are stable since higher-dimensional Schwarzschild 
black hoes are stable [291 EH] • Thus, at critical values of angular momenta, the stability changes and 
there should be a static perturbation. The static perturbation indicates existence of a new family of 
solutions which bifurcates from that of the even- dimensional MP black holes. We found instability in 
the gravitational scalar modes with several k satisfying k > 2. These perturbations break all of the 
symmetry in CP n_1 [5] and, thus, the branched solutions will have only U(l) x R t isometry. To confirm 
this picture by constructing non-linear solutions of the full geometry would be another direction of the 
future research. 

We also found that the mass of the effective field on AdS2, namely, the eigenvalues of the angular 
part of the perturbations are given by simple rational expressions. This property can be found also in 
the four and five-dimensional extreme black holes and odd- dimensional MP black holes [121 HH1 HSl [20]. 
These results prompt us to conjecture that the eigenvalues for the angular part of the perturbations 
are given by rational numbers for any vacuum near-horizon geometries with vanishing cosmological 
constant, as long as we assume the perturbations to be axisymmetric. Further studies on more general 
near-horizon geometries will be useful to falsify such a conjecture. One possible recipient for such an 
analysis is the near-horizon geometries of the odd-dimensional MP black holes with all but one angular 
momenta are equal, which have a cohomogeneity-1 near-horizon geometries. The method used in this 
paper is directly applicable to these geometries, while it involves fifteen unknown variables in general. 
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As a byproduct of our analysis, we found that the expressions of the eigenvalues are essentially 
governed by the transformation property of the fields; Other than a few exceptions, the eigenvalues of 
the scalar field coincide with a part of those of the electromagnetic scalar modes, and the latter coincide 
with a part of those of the gravitational scalar mode. We also found that the eigenvalues have simple 
organization described by two integers (C,£q). We find similar properties also for the vector modes, 
whose eigenvalues are described by Eq. (TT2"j) . It would be interesting to clarify the mathematical and 
physical background for these properties. It would also be nice if we can find microscopic interpretation 
for them based on the dual CFTs. 
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A Background variables 

In the appendices, we show the details of the method to solve the eigenvalue equations In this 
App. [Al we introduce the basic objects which appear in the perturbation equations (JT])-©. Using 
the formulae in this appendix, we show the explicit forms of the perturbation equations and show the 
analytic solutions for them when they exists in App. |B] In the cases we cannot find analytic solutions, 
we solve equations numerically to obtain the eigenvalues. In App. [Cl we summarize the numerical 
technique we employed. 

Below, we work in the frame adapted to the induced metric on the horizon fl3]) given by 

(e e ) M = Q(d9)„ , (e^ = + A)n , (e Q V = #(e Q )^ , (38) 

where e a are real orthonormal frame for CP™" 1 . We show the tetrad components of each tensor quantity 
below unless otherwise noted. The spin connections (c^aft)^ = (e a ) l 'V At (efe) l/ for this frame are given by 

(log$)' (logtf)' $ $ 1 - 

ue<f> — q — e^, ue a — — — e a , — ^J^e^, u a p — — -^Jape^-^ — u a p. (39) 

Assuming the axisymmetry in our case, we have 

k = ft$e^, km 1 = 0, k v k v = fi 2 $ 2 , 

/2$ , $ 2 \ 26 , (40) 

dk = Q I — (log $) e e Ae (j} + -^J af) e a AeJ, dL = — - - (log 6) e e . 
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Using IZab = dujab + u ac A oo cb and de a + oo a b A e& = 0, we find the curvature 2-forms are given by 

n °* = -2^2 ( 2 ( lo s $ )" + ( lo § $ )' ( lo s |) ) e ° A e * - ( log 1) X/3ea A 6/3 ' (41) 

K»a = -zL f 2 ( lQ g - ( lo § 6°S -11 e e A e Q - flog |1 A e p , (42) 



26 2 V V « 

$ / $y /$2 (io g $)'(iog*)'\ 



' 2 1 rh2 



$ / $y (log*)' i - $ 2 / v 

= _ \J/2Q I g ^ J l ^ Q ' 9 661 A e * @2 e ° A e /3 + ^2^ a ' 3 ~~ g>4 WapS/S + Ja [7 1 J^S | <5] J ^Aej , 

(44) 

where the Riemann and Ricci tensors of CP n_1 are given by 

Ra^S = Sa-ySpS ~ ^aS^-y + Ja-yJp& ~ JasJp-y + ^JafiJ-fi , Ra/3 = 2n5 a/ 3 . (45) 



B Perturbation equations and eigenvalues 

We sketch the procedure to construct the eigenvalue equations for X s in this appendix. The proce- 
dure is parallel to that for the odd- dimensional MP black holes with equal angular momenta p2]. We 
show the analytic solutions and eigenvalues for them when they exists, and numerical ones otherwise. 
We defer a part of the derivations of the eigenvalue equations and most of the technical details to 
solve them to App. O The calculations to obtain the tetrad components of the eigenvalue equations 
mentioned below are conducted partially with the aid of the computer algebra system Cadabra [3Tj 132] . 



B.l Scalar field perturbations 

Since our interest is in the axisymmetric perturbations, we take the separation ansatz as 

^ = Xo(t,r)f(9)Y(x a ), (46) 
where Y is the scalar harmonics on CP n_1 defined with respect to the covariant derivative on CP n ~\ 

(D 2 + A;?)Y = 0, \ s K =AK(K + n-l). (k = 0,1,...) (47) 
Decomposing Eq. (j7|) into 6 and CP n_1 parts and using Eq. (]47[) . we find the eigenvalue equation to be 

A /Y = 0<°>/Y = {-<9 2 /Y - (log^a 2 ^))' d e fY - ^D 2 fY + 6 2 M 2 /y} 

- {-r -<*-u - • /' + m + /} ». m 

where /' = df/d9. Note that the parameter a in the metric only appears as a multiplicative factor 
of the mass term in Eq. (1481) . This property originates from the fact that Ao is dimensionless and thus 
only the dimensionless combination of the quantities, aM, may appear on the right-hand side. 
When the scalar field is massless (M = 0), Eq. (1481) has an exact solution given by 

f{9) = dP^cosB) + C^sin^-^costf), (49) 
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where P„{x) and Q£(x) are the associated Legendre functions and 



/i=2K + n — l, u= 2 

The regularity of the solution at 9 = requires C2 to be zero and v to be an integer equal to or larger 
than /i. This requirement gives the eigenvalue Ao shown in Eq. (122]) . 

When M ^ 0, we need to calculate Ao numerically. From Eq. (I48j) . we find that a regular solution 
near 9 = behaves as / ~ 9 2k . Behavior near 9 = ir is the same after 9 t— > rr — 9 due to the symmetry 
about 9 = tt/2 of the near-horizon geometry ([I]). Then, defining a new variable / by / = /sin 2K # 
and imposing the Neumann boundary condition at 9 = and tx to /, we may obtain regular numerical 
solutions by the relaxation method. We defer the further details of the numerical method to App. O 

We show the numerical value of Ao for M = and k = along with the exact value given by 
Eq. (j22j) in Table [21 which shows the relative error of our numerical results from the exact results is 
maintained to be < (9(10 -5 ). We also show the result for aM = 1 and k = in Table [3J As we can 
see from these tables, Ao tends to increase as M increases. As far as we have checked numerically, this 
property holds for any M, k and n. If it holds in general, there is no instability for any M 2 > 0. 



-1 + t/-(2n-l){2(2« + n- l) 2 - 2 + n (4A + 3)} . (50) 
n 



Table 2: Numerical and exact values of Aq for M = and k = 0. 



(a) Numerical (b) Exact (Eq. (1221) ) 



mode# 


n - 2 



1 


2 


3 


4 




mode# 


n - 2 



1 


2 


3 


4 


1 


0.00001 


0.00001 


0.00001 


0.00001 


0.00001 




1 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


2 


1.33330 


1.20001 


1.14287 


1.11112 


1.09091 




2 


1.33333 


1.20000 


1.14286 


1.11111 


1.09091 


3 


3.33334 


2.80000 


2.57143 


2.44445 


2.36364 




3 


3.33333 


2.80000 


2.57143 


2.44444 


2.36364 


4 


5.99997 


4.80001 


4.28573 


4.00001 


3.81818 




4 


6.00000 


4.80000 


4.28571 


4.00000 


3.81818 


5 


9.33333 


7.20001 


6.28572 


5.77778 


5.45455 




5 


9.33333 


7.20000 


6.28571 


5.77778 


5.45455 


6 


13.33330 


10.00001 


8.57144 


7.77779 


7.27273 




6 


13.33333 


10.00000 


8.57143 


7.77778 


7.27273 



Tab 


e 3: Numerical value 


f A for aM = 1 and k = 0. 


mode# 


n - 2 



1 


2 


3 


4 


1 


0.176311 


0.068203 


0.036151 


0.022391 


0.015231 


2 


1.586020 


1.306218 


1.202036 


1.149000 


1.117304 


3 


3.600436 


2.920479 


2.641832 


2.491085 


2.396961 


4 


6.271692 


4.927228 


4.362393 


4.052046 


3.856095 


5 


9.607142 


7.330994 


6.366260 


5.833399 


5.495654 


6 


13.608308 


10.133272 


8.654517 


7.835879 


7.316149 



B.2 Electromagnetic perturbations 

Next, we discuss the perturbations of the Maxwell field on the near-horizon geometries. For the 
analysis, it is useful to decompose the perturbation variables Y" M into the scalar and vector modes 
according to the transformation properties with respect to CP™" 1 . We show the derivations of the 
eigenvalue equations and results mode by mode below. 
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B.2.1 Electromagnetic vector modes 

The electromagnetic vector modes are characterized by 

y 4 = 0, D ±a Y a = 0, (i = 6,(f>) (51) 

where D ± is the covariant derivative on CP n_1 projected to the ^pi eigenspaces of the complex structure 
for CP n_1 , J = \j a e a efi, defined as 

Dt = V^Dp, V% = \ (g a p ± iJap) . (52) 

We may parametrize this component of the perturbations as Y a = f(9)Y a , where Y a is the divergence- 
free vector harmonic satisfying 

D ±a Y a = 0, Y± = -^=Df a Yp } , (D 2 + X v K )Y a = 0, X y K = 4k(k + 2) + 2n (2k + 3) . (53) 

For these modes, we find i components of Eq. (JHj) becomes trivial, and CP n_1 components reduce to an 
equation of f(9) given by 

n(9n m m\-f"((j\ i {2n - l){(2n + l)cos 2 g- 1} 
- n(2n - l)X 1 f(9) - f (9) + — {( ^ _ ^ _ - i} / (61) 



+ 



(2n- l) 2 cos 2 9 + An - 1 (2k + 3)(2n + 2k + 1) 



{(2n — 1) cos 2 6> + l} 2 sin 2 # 



| 2(2n-l)((n + 2)K + n + K 2 ) l (g4) 
n J 

This equation has an analytic solution given by 

/(0) = dPftcoaO) + ^(sinfl^Q^cosfl), (55) 

where 



/x = 2k + n + 2, ^ = 2 



-1 + \l - {4n 3 + 4n 2 (4K + 2A X + 5) + ti(8k(2k + 3) - 4Ai - 7) - 8k(k + 2)} 
n 



(56) 

For the solution to be regular at 9 = 0, C*2 should be set to zero and v should be an integer equal to 
or larger than \x. This requirement fixes Ai as given in Eq. ( 1231) . 

B.2.2 Electromagnetic scalar modes 

Next, we discuss the electromagnetic scalar modes for which % components of the perturbations are 
turned on. We may expand the scalar mode perturbations as 

Y = f t (9)Y, Y a = g + (9)Yt + g~(9)Y~, (57) 

where are the scalar-derived vector eigenfunctions defined by 
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We further assume that are real-valued, and introduce the real quantities by 

g ± = g R ±ig I , Y± = Y* ± zY J a . (59) 

As a result, we find coupled ODEs for (f e , f^, g R , g 1 ) as the eigenvalue equations when k > 0. For 
k = 0, Y becomes a constant and the variables reduce to (fg, f^). We solve those equations using the 
numerical technique demonstrated in App. O 

In Table HJ we show an example of the numerically-obtained eigenvalues Ai. As far as we have 
examined, our numerical result is consistent with a hypothesis that the eigenvalues are expressed in 
general by Eqs. (Ml, (123]) and (EE). 



Table 4: Numerical value of Ai of the electromagnetic scalar modes for k = 4 and n — 2, . . . , 6. 



K = 4 


n-2 



1 


2 


3 


4 


mode# 


1 


4.000000 


1.600000 


0.857143 


0.533333 


0.363636 


2 


9.333331 


4.799999 


3.142856 


2.311111 


1.818181 


3 


9.999997 


5.599998 


3.999999 


3.199999 


2.727272 


4 


9.999997 


5.599998 


3.999999 


3.199999 


2.727272 


5 


15.999992 


9.199996 


6.571425 


5.199997 


4.363634 


6 


16.666658 


9.999995 


7.428567 


6.088885 


5.272724 


7 


16.666656 


9.999994 


7.428566 


6.088885 


5.272723 


8 


16.666658 


9.999995 


7.428568 


6.088885 


5.272724 


9 


23.333315 


13.999989 


10.285706 


8.311104 


7.090903 


10 


23.999980 


14.799987 


11.142846 


9.199990 


7.999991 


11 


23.999983 


14.799990 


11.142849 


9.199993 


7.999994 


12 


23.999981 


14.799988 


11.142847 


9.199991 


7.999992 



B.3 Gravitational perturbations 

Finally, we study the gravitational perturbations. The procedure is parallel to that for the electro- 
magnetic perturbations, while the gravitational perturbations involve ten unknown variables in general. 
We show the procedures and results for each of tensor, vector and scalar modes below. 

B.3.1 Gravitational tensor modes 

We start from the simplest components, which are the gravitational tensor modes defined by 

5^ = 1^ = 0, g af> Y a p = Q, D ±a Y aW =0. (60) 

Under these assumptions, only the CP n_1 components of Eq. ([9]) remain nontrivial. We introduce a 
decomposition given by 

Y a(3 = f(9)Y al3 , (61) 

and we also decompose Y a p into hermitian and anti-hermitian parts by {JYJ') a!3 = crY a/ 3, where 
a = — 1 (+1) for the hermitian (anti-hermitian) modes. 
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To deal with the D 2 Y a p terms in the equations, we need to take Y a p to be an eigenstate of the 
generalized Lichnerowicz operator on CP n ~ , that is, 



\lY aP = (A^Y) q/3 = -D 2 Y a p - 2R a7/3S Y' yS + 4nY a/3 = -D 2 Y a p + 2{2n + l)Y a p - 3 (JYJ) af3 . (62) 

For n = 2, the tensor harmonics do not exist on CP 1 = S 2 . For n > 3, the harmonics exist and the 
eigenvalues of are given by [TU] 



Xl = 4k (k + n - 1) + 4 (n - 1 + a) . 



(63) 



This allows us to rewrite D 2 Y a p in terms of Y a p. Using the above equations, the eigenvalue equation 
flU is rewritten as 



(2n-l)A 2 /(0) = /"(0) + (2n-l) 



2n 



cot9 



2 sin 29 



(2n - 1) cos 2 9 + 1 



f'(0) 



+ K-1 - 



2k(h + k-1] 



(2n - 1) cos 2 9 + 1 ' " V sin 2 # 
-k 2 + K — a + 1 



+ 



n 



k(2k - 3) + 2a - 3 



This equation has an exact solution given by 

f(0) = dP^cosfl) + C 2 (sm0) 1 - n gC(cosi 

where 



/(*)■ (64) 



(65) 



i / I2n — \ 

H = 2k + n - 1, i/ = - ( -1 + ^/ — — {2 {2k + n - l) 2 - 10 + n (4A 2 + 3) + 8a} ) . (66) 



n 



Setting C*2 = and z/ > \i to be an integer for the sake of the regularity at 9 = 0, we find A 2 to be 
given by Eq. (TIT]) . 

B.3.2 Gravitational vector modes 

The gravitational vector modes are composed of divergence free vectors Y ia and traceless Y a p sat- 
isfying 

Yy = 0, D ±a Y ta = 0, Y a a = 0. (67) 
We expand the perturbations as 



Y ia = 9i (e)Y a , Y a ? = h + (9)Y^ + h-{9)Y~ p ee Y+ + 3^, 



(6* 



where Y a is a divergence- free vector harmonics defined by Eq. fl53l) . We further decompose Y a into the 
eigenvectors of J a p by J^Yp = —ieY a with e — ±1. 

Since Y a are complex, the variables (g^, (/i + + h~)/2, (h + — h~)/2i) are not real-valued. To rewrite 
the variables in terms of real quantities, we further decompose the variables as 



Y ia = g+(0)Y+ + g-(9)Y~, Y afi = h+ + (9)Y++ + h- + (9)Y~; + h+~(9)Y+- + h-(9)Y~- (69) 
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where Y^ are the vector harmonics corresponding to e 

-1 



±1 (that is, J a 



7 + 



=RY±) and 



±1±2 



Y^ 1 



Then, replacing the quantities by 

9t = 9? ± wL 
Yf = Y R ± iYj, 



,±1± 2 



11 a/3 



(70) 



(h RR ± 2 ih RI ) ± 1 i(h IR ± 2 ih H ), 

(Y«?± a iYS) ±ii(Y3± 2 iY5), 



(71) 



we obtain real- valued eigenvalue equations in terms of (g R , g\, g R , g^, h RR , h RI , /i 7ii , /i JJ ). 

Due to the self-adjointness of the operators with respect to suitable inner products, A2 is 
guaranteed to be real numbers [12J. This fact implies that the e = ±1 modes share the same eigenvalue 
A 2 . 

We show the numerically-obtained eigenvalues in Table [51 Note that this result contains the eigen- 
values for both e = 1 and e = — 1 modes, which are equal to each other. This numerical result 
suggest that there are four species of the eigenvalues, other than the multiplicity coming from e = ±1 
modes, and the expression of the eigenvalues for general n and k is given by Eqs. (IT2"|) and (fT3]) . These 
eigenvalues are all positive, and hence no instability is implied. 



Table 5: A2 of gravitational vector modes for n = 5 and k = 0, . . . , 3. 



n = 5 


K 








mode# 





1 


2 


3 


1 


0.4444 


0.5778 


0.8000 


1.1111 


2 


0.4444 


0.5778 


0.8000 


1.1111 


3 


0.8889 


1.4667 


2.1333 


2.8889 


4 


0.8889 


1.4667 


2.1333 


2.8889 


5 


2.0000 


2.5778 


3.2444 


4.0000 


6 


2.0000 


2.5778 


3.2444 


4.0000 


7 


2.0000 


2.5778 


3.2444 


4.0000 


8 


2.0000 


2.5778 


3.2444 


4.0000 


9 


2.6667 


3.6889 


4.8000 


6.0000 


10 


2.6667 


3.6889 


4.8000 


6.0000 


11 


3.7778 


4.8000 


5.9111 


7.1111 


12 


3.7778 


4.8000 


5.9111 


7.1111 



n = 5 


K 








mode# 





1 


2 


3 


13 


3.7778 


4.8000 


5.9111 


7.1111 


14 


3.7778 


4.8000 


5.9111 


7.1111 


15 


3.7778 


4.8000 


5.9111 


7.1111 


16 


3.7778 


4.8000 


5.9111 


7.1111 


17 


4.6667 


6.1333 


7.6889 


9.3333 


18 


4.6667 


6.1333 


7.6889 


9.3333 


19 


5.7778 


7.2444 


8.8000 


10.4444 


20 


5.7778 


7.2444 


8.8000 


10.4444 


21 


5.7778 


7.2444 


8.8000 


10.4444 


22 


5.7778 


7.2444 


8.8000 


10.4444 


23 


5.7778 


7.2444 


8.8000 


10.4444 


24 


5.7778 


7.2444 


8.8000 


10.4444 



B.3.3 Gravitational scalar modes 

Finally, we comment on the gravitational scalar modes, for which the analysis procedure is similar 
to the previous examples while it is more involved. 

We expand the metric perturbations using the scalar harmonics Y as 



Y ia = g+(9)Y+ + g7(9)Y-, 
1 



(72) 



Yap = — 7p {h ++ mU + h-(e)Y-p + h + -(6)Yip) - ——- (f ee (6) + / # (0))<5 a/3 Y, 



where Y^ and Y Q ^ are scalar-derived vector /tensor eigenfunctions defined by Eq. ( 1581) . To construct 
the real-valued equations, we redefine the variables and mode functions as 

g± = g*±ig I i , h ±± = h R ± i h I , Y± = Y«±zY^, Y±± = Y*±zY^. (73) 

As a result, we obtain coupled ODEs of (fgg, fg^, f^, g R , gjj, g R , gl, h R , h 1 , h + ~) in a general case (k > 1 
and n > 2) as the eigenvalue equations to solve. Depending on the values of k and n, some of the 
unknown variables drop out as follows (see also Ref. [T2]): 

• k = 0: Y becomes a constant and there is no scalar-derived vector nor tensor in this case. We 
have only (/y) as the unknown variables. 

• K = l& n = 2: All of Y^ and Y^~ vanish in this case, and we have (/y, gf' 1 ) as the unknown 
variables. 

• K=l&n>2: Y^ vanish in this case, and we have (fijyg^ 1 , h + ~) as the unknown variables. 

• k > 1 & n = 2: Y^~ vanishes in this case. As a result, we have (/y, gf 1 ' 7 , /i ±=t ) as the unknown 
variables. 

The numerically-obtained eigenvalues for k > 1 and n > 2 suggest that they are described by 
Eqs. (I14p and f lT5|) . We show some examples of the numerically obtained eigenvalues A2 in Table [6j For 
the special cases mentioned above, we find that the eigenvalues are given by Eq. (114]) in any cases, and 
only some of the integer sets (C, £q) are modified. We list the integer sets for each case in Eqs. (|16p - 
(]T9"j) . In short, some of the eigenvalues listed in Eq. (|15p vanishes, and also £0 is shifted in some of the 
remaining ones. 



C Technical details 

We summarize the numerical technique to solve the eigenvalue equations in this appendix. This 
technique is essentially the same as that used in Ref. [18], which studied the conformal weights for the 
general near-horizon geometries in five dimensions. 

C.l Perturbation equations and boundary conditions for electromagnetic 
scalar modes 

We explain the numerical technique to find the eigenvalues taking the electromagnetic scalar modes 
for example. For the electromagnetic scalar modes defined by Eq. (15 7p . we obtain a system of coupled 
ODEs for the unknown variables (fg, f$, g R , g 1 ) if k > 0, where these variables are defined by Eqs. (157]) 
and ( )59|) . If k = 0, the scalar harmonics Y becomes a constant, and the unknown variables reduce to 
(fe, f(j>)- 

Below, we take the k > case as an example, and explain the boundary conditions to be imposed 
on the unknown variables. The resultant eigenvalue equations may be expressed as 

(d 2 + A(d)d e + B{d)) v = -(2n - 1)A^, (74) 

where v = [fg, fag*^ 1 ) , and A{6) and B(9) are 4x4 coefficient matrices which behave near 8 = 

as 

A= 2 -^-±E + O(9 ), B = ±B_ 2 + 0(9^), (75) 
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Table 6: A2 of gravitational scalar modes for n = 2, 5 and n = 0, . . . , 4. 



77, — Z 


K 












n 
u 


1 
1 


9 


Q 
O 


A 


1 


3.333 


1.333 


-0.667 


-0.667 


-0.000 


2 


4.000 


2.000 


0.667 


2.000 


4.000 


3 


5.333 


2.000 


2.000 


3.333 


5.333 


4 


6.000 


3.333 


2.000 


3.333 


5.333 


5 


6.667 


4.000 


3.333 


6.000 


9.333 


6 


8.000 


4.667 


4.000 


6.667 


10.000 


7 


9.333 


4.667 


4.000 


6.667 


10.000 


8 


10.000 


4.667 


5.333 


8.000 


11.333 


9 


11.333 


6.000 


5.333 


8.000 


11.333 


10 


13.333 


6.000 


5.333 


8.000 


11.333 


11 


14.000 


7.333 


7.333 


11.333 


16.000 


12 


15.333 


8.000 


8.000 


12.000 


16.667 


13 


18.000 


8.000 


8.000 


12.000 


16.667 


14 


18.667 


8.000 


8.000 


12.000 


16.667 


15 


20.000 


9.333 


9.333 


13.333 


18.000 


16 


23.333 


9.333 


9.333 


13.333 


18.000 


17 


24.000 


9.333 


9.333 


13.333 


18.000 


18 


25.333 


11.333 


9.333 


13.333 


18.000 


19 


29.333 


12.000 


12.000 


17.333 


23.333 


20 


30.000 


12.000 


12.667 


18.000 


24.000 



71 — 5 


K 












n 

u 


1 
1 


9 


Q 
O 


A 


1 


2.444 


0.444 


-0.356 


-0.622 


-0.800 


2 


3.333 


1.333 


0.089 


0.267 


0.533 


3 


4.000 


1.333 


0.978 


1.378 


1.644 


4 


4.444 


2.000 


1.200 


1.378 


1.644 


5 


4.889 


2.444 


1.200 


1.600 


2.311 


6 


5.778 


2.444 


1.867 


2.489 


3.200 


7 


6.000 


2.889 


1.867 


2.489 


3.200 


8 


6.667 


2.889 


2.978 


3.600 


4.311 


9 


7.778 


2.889 


2.978 


3.600 


4.311 


10 


7.778 


3.778 


2.978 


3.600 


4.311 


11 


8.667 


4.000 


2.978 


3.600 


4.311 


12 


9.778 


4.000 


2.978 


4.044 


5.200 


13 


10.000 


4.000 


3.867 


4.933 


6.089 


14 


10.889 


4.667 


3.867 


4.933 


6.089 


15 


12.000 


4.667 


3.867 


4.933 


6.089 


16 


12.444 


4.667 


4.978 


6.044 


7.200 


17 


13.333 


5.778 


4.978 


6.044 


7.200 


18 


14.444 


5.778 


4.978 


6.044 


7.200 


19 


15.111 


5.778 


4.978 


6.044 


7.200 


20 


16.000 


5.778 


4.978 


6.044 


7.200 



where E is the unit matrix. To find out the fall-off behavior of the solution at 9 = 0, we need to 
diagonalize the coefficient matrix B_ 2 by introducing new variables v by v = Su, where S is the matrix 
composed of the eigenvectors of B_ 2 . The equations in terms of these new variables are 

(d 2 9 + A(9)d e + B{&)\ v = -(2n- 1)A^, (76) 
where A = + O(6 ) and B behaves for 9 -> as 

B = ^B- 2 + 0(6- 1 ), B_ 2 = Diag [e_, e_, e+, e+] , e± = - (2« ± 1) {2(n - 1) + 2« ± 1} . (77) 
Solving Eq. (J75J) near 6 = 0, we find the asymptotic behavior of the solution near 9 = to be given by 

* m \ 2k ~ 1 ( i = 1 ' 2 ) (7^ 

Vi oc 9 P \ Pi = < , 78 

F (2k + 1 (i = 3,4) V ; 

where we chose the decaying solutions from the two linearly-independent solutions so that the solutions 
are regular at 9 = 0. Based on this observation, we introduce new variables v by v = Tlv, where II is 
the diagonal matrix whose components are given by Iljj = sin Pi 9. For these new variables v, we find 
the eigenvalue equations become 

(d 2 e + A(9)d e + B{6)\ v = -(2n - l)Xiv, (79) 

where A = 0(9~ l ) and B = O(9 ). This equation has a regular singularity at 9 = 0, and we need to 
impose the Neumann boundary condition there, namely, v = vq + 0(9 2 ), to maintain the regularity of 
the solutions. 
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Since the perturbation equations (jHJ) become symmetric with respect to 9 = n /2 for the near-horizon 
metric (pQ), we may solve Eq. ( !79|) by imposing the Neumann boundary condition at both 6 = and 
9 = it. Alternatively, we separate the variables into the even/odd modes with respect to 9 = n/2 as 
follows. In Eq. (!74|) . we find the parity with respect to 9 = n/2 of the coefficient matrix components 

1 i=2,3,4 



are given as 

A = i=1 { 

i=2,3,4{ 



(odd) 


(even) 


(even) 


(odd) 



B 



(even) 


(odd) 


(odd) 


(even) 



30) 



This property originates from that the original equations are symmetric with respect to 9 = n/2, and 
also that the variable V\ = fg is an odd function with respect to 9 = n/2, since it is proportional to 
the eg component of the perturbation variables, while 1*2,3,4 are not. This property implies that an 
eigenvector take a form either of 





where v\ and vf 1 are odd and even functions with respect to 9 = n/2, respectively. Since both 
A and B are regular at 9 = n/2, we may impose the Dirichlet boundary condition on v\ to set it 
zero and the Neumann boundary condition on ^2,3,4 to obtain solutions corresponding to the left one 
in Eq. (18 ip . Solutions corresponding to the right one in Eq. (jSTj) can be obtained by switching the 
boundary condition types. By decomposing into the even/odd modes, we may reduce the calculation 
region < 9 < n to the half, and also reduce the number of modes needs to be subtracted at each 
time step in the relaxation method, which will be introduced in the next section, to the half. It helps 
increasing the precision of the numerical results in a shorter calculation time. 

In practice, the even/odd decomposition is more simply implemented by introducing new variables 
v by v = Tv, where T is a diagonal matrix whose components are equal tqj cos(^(l — COS6 1 )) for i = 1 
(i = 2, 3, 4) and 1 for i = 2, 3, 4 (i = 1) if we are to find the solutions of the left (right) type in Eq. (|8ip . 
In terms of v, the eigenvalue equations become 

(d 2 e + A(8)d e + B{6)) i = -(2n - l)Xiv, (82) 

where A = — 9)~ x ) and B{9) = C((f — 9)°), and we may impose the Neumann boundary 

condition on all the components of v to find the solutions of the left (right) type in Eq. ( IHTi) . Following 
the procedures in the previous paragraphs after introduction of T, we define and use the variables given 
by v = (TSn) -1 ^, and solve Eq. (179|) imposing the Neumann boundary condition at 9 = and n/2. 

For k = 0, the unknown variables reduce to v = (fg, f^) as we mentioned previously. We find B_ 2 = 
Diag[— (2n — 1), — (2n — 1)] without the diagonalization, and we find that the regular solutions behave 
as Vi oc 9°. Thus, we may solve the eigenvalue equations using the original variables Vi, or alternatively 
introducing only the even/odd decomposition defined by T = Diag(w, 1) or T = Diag(l,u>) with 
w = cos( : |(l — cos 9)) . 



4 The function w — cos(?(l — cos6>)) is chosen so that w and w/(6 — n/2) are sufficiently close to constants near 
= and 6 = n/2, respectively. For example, simpler choices such as w = cos# or w = cos(26 2 /ir) give rise to 0(9~ 2 ) 
or 0[(9 — f )~ : ) terms in the components of the matrix B in some cases. Such terms in B become obstacles for our 
scheme. 
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C.2 Numerical implementations 

We solve the eigenvalue equation (179|) numerically as follows. We employ the relaxation method, 
which was used also in Ref. [18] to treat similar problems. 
Firstly, we introduce a diffusion equation given by 

dj( r , 9) = M(9)v, M(6) = dj + A{6)d e + B(9), (83) 

whose solution behaves as 

£( r ,0) = e rM d init -> e T H, (84) 

where A* is the largest eigenvalue of M and v* is the corresponding eigenvector. It implies that we may 
obtain the largest eigenvalue by following the time evolution described by Eq. (183]) for a sufficiently 
long time. The smaller eigenvalues can be obtained successively by the same calculations if we project 
out the eigenvectors for larger eigenvalues at each step of the time evolution. 

To make the equation amenable to numerics, we need to discretize the diffusion equation ( 18"3"|) on 
the grids given by r = r n = nAr (n = 0, 1, . . .) and 9 = 9 { = iA9 (i = 0, 1, . . . , N, Ad = n/2N). We 
do it implicitly as, for 1 < % < N — 1, 

~~^T~ = A^ + 2A9 + B ^ ■ (85) 

An advantage of the implicit scheme is that the time evolution is stable for any At, and we may take 
it large to speed up the numerical calculations. Eq. ( [85]) may be transformed into 



V?' 1 



= a^_i + ht? + Ci^ +1 . (86) 

We need to impose the Neumann boundary condition at 9 = and 9 = tt/2. At 9 = 0, it implies 
v ~ v + \v'q9 2 , and Eq. fl83|) at 6 1 = may be expressed as 

d T i> = (e + Ao) d 2 e £> + B £>, (87) 
where A Q = lim^o This equation is discretized to give an equation for z = as 



Ar A9 



Following the same procedure for 9 = tc/2 , we have the equation at i = N as 
2Ar / ? \ // 2Ar\ „ 2Ar 



,,n-l 
N 



A9 2 

where = lim fl _>2 (0 - §)i(0). 



+ Iat) + ((l + + ^T^v - AtB n ^) d n N = aJ n N _, + 6^, (89) 
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The above equations are summarized as follows: 



-n-l 



-n-l 
'N-l 
-n-l 



( b c 
ai bi ci 



\ 



V 



a N-l &ZV-1 OV-l 

Q>N bft J 



v" 
' 



U N-1 

\v n N ) 



(90) 



This equation is readily solved by the forward and backward substitution to obtain v n from v n ~ l . To 
obtain the full spectrum, we need to conduct two calculations for the two T matrices, each of which 
corresponds to the each mode in Eq. ( IHTj) . 

Applying this method to Eq. (!79|) for the electromagnetic scalar modes, we obtain the results 
summarized in Table H] and also as Eq. ( 124")) supplemented with Eq. (125|) . Similarly, we may obtain 
the results for the k = case, and the results are summarized as Eq. (1261) . In this calculation, we 
took iV = 4000 and At = 10 -2 , and truncated the time evolution once the condition |A _1 ^| < 10 -8 is 
achieved. After the truncation, we progressed the time evolution one more step using smaller At = 10 -8 
to obtain the correct value of Ai. 

C.3 Gravitational vector modes 

The gravitational vector modes have the following properties. The unknown variables are 







R „I „R „I faRR frRI uIR uII\ T 



(91) 



whose definitions are found in Eq. (I7TT) . Following the procedure of Sec. IC.ll in this case, we find the 
diagonalized coefficient matrix B_ 2 is given by 



Diag [e_, e_, e_, e_, e+, e+, e+, e+] , 



-4(«+2)(n+«), 



-4(K+2)(n+K+l). (92) 



Solving the equation corresponding to Eq. (1761) . we find the asymptotic behaviors for 6 — > of the 
decaying solutions to be 



6 p -,6 p -,6 p +,9 p+ ,9 p+ i6 p+ ] 



p_ = 2(« + l), P+ = 2(k + 2). 



(93) 



Then, we define and use new variables v = U~ 1; v, where II is a diagonal matrix whose components are 
given by Iljj = sin p ~ 9 for % = 1, . . . , 4 and Ua = sin p+ 9 for i = 5, . . . , 8. Following the procedure in 
Sees. IC.ll and IC.2t we obtain the eigenvalues summarized in Tableland also as Eqs. f|T2l) and (TT3l) . 



C.4 Gravitational scalar modes 

The gravitational scalar modes have the following properties. For a general case with k > 1 and 
n > 2, the unknown variables are given by 



v = (fee, fe<t>, fn>, 9e, 9e> 9*, 9d>, h R , h 1 ', h + ) 



T 



whose definitions are found in Eq. (173]) . B_ 2 for this case becomes 

iL 2 = Diag [e_, e_, e_, e , e , e , e , e+, e+, e+] , e s = -4 (k + 5) (k + n - 1 + 5) 



(94) 



(95) 
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which give rise to the asymptotic behavior described as 

v~ (6 p -,6 p -,6 p -,6 P0 ,6 p0 ,6 P0 ,6 P0 ,6 p+ ,6 p +,6 p+ f , Ps = 2{k + 5). (96) 

We define and use the matrix IT using p$ in this equation. The numerical results of the eigenvalue A2 
are summarized as Eq. (124)) with Eq. (125]) . 

Below, we comment on the special cases for the gravitational scalar modes. 

• k = 0: 

Y = constant in this case, and the unknown variables are given by v = {foe, fg^, f^). Without 
introducing the diagonalization, we find _B_2 = Diag[— 4n, — An, — An], from which we find that 
the regular solutions behave as oc 9 2 . We define and use E = E and IT such that 11^ = sin 2 9. 

• K=l&n = 2: 

All of and Y+^ vanish, and as a result we have only (f ee , f d(t> , f ' H , gf, g\, gf, g 1 ^) as the free 
variables in this case. It follows _B_2 = Diag[0, —8, —8, —8, —24, —24, —24], from which we find 
Vi oc 9 Pi with pi = 0, P2,3,4 = 2 and 255,6,7 = 4. 

• k = 1 k n > 2: 

Y^ vanishes in this case, and we have v = {fee, fe<t>, f^i 9^ 9di 9$i 9\i h + ~) T as the free variables. 
We find 5_ 2 = Diag[0, — An, —An, —An, —8{n + 1), — 8{n + 1), — 8(n + 1)], which results in Vi oc 9 Pl 
with pt = 0, p 2 ,3,4 = 2 and p 5A7 = A. 

• k> 1 & n = 2: 

Y^T vanishes in this case, and as a result we have {fee, fee/,, /</,</>, gf, gj), g^, gff,,h ++ , h ) as the 
free variables. We find B_ 2 = Diag[e_, e_, e_, e , e , e , e+, e+, e+], where e$ are those in Eq. (155]) 
with n = 2 plugged in. The asymptotic behavior of the regular solutions in this case is given by 
Vi oc 9 Pi with pi,2,3 = 2(k - 1), f> 4 ,5,6 = 2k and p 7 , 8 , 9 = 2(ac + 1). 

The numerical results of the eigenvalues A2 for these cases are summarized by Eqs. (fIB"]) - (lT9~]) . 



References 

[1] R. Emparan, "Black holes galore in D > 4," Fortsch. Phys. 56, 723 (2008). 

[2] R. Emparan and H. S. Reall, "Black Holes in Higher Dimensions," Living Rev. Rel. 11, 6 (2008) 
|arXiv:0801.347T1 [hep-th]]. 

[3] R. C. Myers and M. J. Perry, "Black Holes in Higher Dimensional Space-Times," Annals Phys. 
172, 304 (1986). 

[4] R. Emparan and R. C. Myers, "Instability of ultra-spinning black holes," JHEP 0309, 025 (2003) 
|hep-th/0308056] . 

[5] O. J. C. Dias, P. Figueras, R. Monteiro, H. S. Reall and J. E. Santos, "An instability of higher- 
dimensional rotating black holes," JHEP 1005, 076 (2010) farXiv: 1001 .45271 [hep-th]]. 



24 



H. Kodama, R. A. Konoplya and A. Zhidenko, "Gravitational stability of simply rotating Myers- 
Perry black holes: Tensorial perturbations," Phys. Rev. D 81, 044007 (2010) jarXiv: 0904 .21541 
[gr-qc]]. 

O. J. C. Dias, P. Figueras, R. Monteiro, J. E. Santos and R. Emparan, "Instability and new phases 
of higher-dimensional rotating black holes," Phys. Rev. D 80, 111701 (2009) jarXiv:0907.2248l [hep- 
th]]. 

M. Shibata and H. Yoshino, "Bar-mode instability of rapidly spinning black hole in higher 
dimensions: Numerical simulation in general relativity," Phys. Rev. D 81, 104035 (2010) 
jarXiv: 1004.49701 [gr-qc]]. 

O. J. C. Dias, P. Figueras, R. Monteiro and J. E. Santos, "Ultraspinning instability of rotating 
black holes," Phys. Rev. D 82, 104025 (2010) [arXiv: 1006. 19041 [hep-th]]. 

H. K. Kunduri, J. Lucietti and H. S. Reall, "Gravitational perturbations of higher dimensional 
rotating black holes: Tensor perturbations," Phys. Rev. D 74, 084021 (2006) |hep-th/0606~0~76~] . 

K. Murata and J. Soda, "Stability of Five-dimensional Myers-Perry Black Holes with Equal An- 
gular Momenta," Prog. Theor. Phys. 120, 561 (2008) iarXiv:0803.137Tl [hep-th]]. 

M. Durkee and H. S. Reall, "Perturbations of near-horizon geometries and instabilities of Myers- 
Perry black holes," Phys. Rev. D 83, 104044 (2011) [ arXiv: 1012.48051 [hep-th]]. 

O. J. C. Dias, R. Monteiro and J. E. Santos, "Ultraspinning instability: the missing link," JHEP 
1108, 139 (2011) jarXiv: 1106.4551 [hep-th]]. 

M. Durkee and H. S. Reall, "Perturbations of higher- dimensional spacetimes," Class. Quant. Grav. 
28, 035011 (2011) (arXiv: 1009.00151 [gr-qc]]. 

M. Durkee, V. Pravda, A. Pravdova and H. S. Reall, "Generalization of the Geroch-Held-Penrose 
formalism to higher dimensions," Class. Quant. Grav. 27, 215010 (2010) [arXiv: 1 002.4826 [gr-qc]]. 

S. A. Teukolsky, "Rotating black holes - separable wave equations for gravitational and electro- 
magnetic perturbations," Phys. Rev. Lett. 29, 1114 (1972). 

M. Guica, T. Hartman, W. Song and A. Strominger, "The Kerr/CFT Correspondence," Phys. 
Rev. D 80, 124008 (2009) jarXiv:0809.4266l [hep-th]]. 

K. Murata, "Conformal weights in the Kerr/CFT correspondence," JHEP 1105, 117 (2011) 
jarXiv: 1103.56351 [hep-th]]. 

O. J. C. Dias, H. S. Reall and J. E. Santos, "Kerr-CFT and gravitational perturbations," JHEP 
0908, 101 (2009) jarXiv:0906.2380l [hep-th]]. 

A. J. Amsel, G. T. Horowitz, D. Marolf and M. M. Roberts, "No Dynamics in the Extremal Kerr 
Throat," JHEP 0909, 044 (2009) jarXiv:0906.2376l [hep-th]]. 

H. K. Kunduri and J. Lucietti, "A classification of near-horizon geometries of extremal vacuum 
black holes," J. Math. Phys. 50, 082502 (2009) jarXiv:0806.2"05Tl [hep-th]]. 



25 



[22] H. K. Kunduri, J. Lucietti and H. S. Reall, "Near-horizon symmetries of extremal black holes," 
Class. Quant. Grav. 24, 4169 (2007) [arXiv:0705.42l4l [hep-th]]. 

[23] S. Hollands and A. Ishibashi, "All vacuum near horizon geometries in arbitrary dimensions," 
Annales Henri Poincare 10, 1537 (2010) |arXiv:0909.3462l [gr-qc]]. 

[24] P. Figueras, H. K. Kunduri, J. Lucietti and M. Rangamani, "Extremal vacuum black holes in 
higher dimensions," Phys. Rev. D 78, 044042 (2008) |arXiv:0803.2998l [hep-th]]. 

[25] S. W. Hawking and H. S. Reall, "Charged and rotating AdS black holes and their CFT duals," 
Phys. Rev. D 61, 024014 (2000) |hep-th/9908109] . 

[26] H. Kodama, "Superradiance and Instability of Black Holes," Prog. Theor. Phys. Suppl. 172, 11 
(2008) |arXiv:0711.4"l8l [hep-th]]. 

[27] P. Figueras, K. Murata and H. S. Reall, "Black hole instabilities and local Penrose inequalities," 
Class. Quant. Grav. 28, 225030 (2011) [ arXiv: 1107. 57851 [gr-qc]]. 

[28] S. Hollands and R. M. Wald, "Stability of Black Holes and Black Branes." ErXiv: 1201 .04631 [gr-qc], 

[29] A. Ishibashi and H. Kodama, "Stability of higher dimensional Schwarzschild black holes," Prog. 
Theor. Phys. 110, 901 (2003) |hep-th/0305l"85] . 

[30] R. A. Konoplya and A. Zhidenko, "Stability of multidimensional black holes: Complete numerical 
analysis," Nucl. Phys. B 777, 182 (2007) |hep-th/070323~I1 [HEP-TH]]. 



[31] K. Peeters, "Introducing Cadabra: A Symbolic computer algebra system for field theory problems, 



hep-th/0701238 



[32] K. Peeters, "A Field-theory motivated approach to symbolic computer algebra," Comput. Phys. 
Commun. 176, 550 (2007) |cs/0608005| [cs.SC]]. 



26 



